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We present new experimental results of low temperature x-ray synchrotron diffraction, neutron 
scattering and very low temperature (mK-range) bulk measurements on the nanotube system 
Na2V307. The crystal structure determined from our data is similar to the previously proposed 
model (P. Millet et al. J. Solid State Chem. 147, 676 (1999)), but also deviates from it in signifi- 
cant details. The structure comprises nanotubes along the c-axis formed by stacking units of two 
V-rings buckled in the afc-plane. The space group is P3 and the composition is nonstoichiometric, 
Na2~a:V307, x=0.17. The thermal evolution of the lattice parameters reveals anisotropic lattice 
compression on cooling. Neutron scattering experiments monitor a very weak magnetic signal at 
energies from -20 to 9 meV. New magnetic susceptibility, specific heat measurements and decay of 
remanent magnetization in the 30 mK - 300 mK range reveal that the previously observed transition 
at «76 mK is spin-glass like with no long-range order. Presented experimental observations do not 
support models of isolated clusters, but are compatible with a model of odd-legged S=l/2 spin 
tubes possibly segmented into fragments with different lengths. 



PACS numbers; 61.66.-f , 73.22.-f, 75.25. +z, 75.50.Lk 

I. INTRODUCTION 

The perspectives to create magnetic and electronic de- 
vices on the atomic scale motivates nowadays research 
on spin clusters, chains, ladders and other geometrical 
spin arrangements. Systems that accommodate such low- 
dimensional spin arrangements in a periodic fashion, i.e. 
on a lattice, are of particular interest. This is because 
their properties arising from individual objects and from 
the whole assemble can be probed by experimental and 
theoretical techniques developed for crystalline solids. 
One of the rare realizations of a crystalline system formed 
by nanotubes with periodic arrangement is Na2V307, a 
compound first synthesized by Millet et alii . The struc- 



ture is formed by VO5 square pyramids connected via 
edges and apices (Fig. [T^). There are four inequivalent 
sodium sites. One Na site is located in the middle of 
the nanotube and the other three around the nanotubes. 
Buckled rings formed by nine V atoms (only three sites 
are inequivalent) have two different orientations forming 
slices A and B (Fig. [TI3). The nanotube is formed by an 
-A-B-A-B- arrangement along the c-axis. 
Extensive experimental investigations employing specific 
heat, ac susceptibility and ^^Na nuclear magnetic reso- 
nance by Gavilano et al^'^ revealed puzzling magnetic 
properties. Those results implied a gradual reduction 
of the number of V"'+ S=l/2 magnetic moments, which 
may be due to the formation of dimers with most of them 
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in a singlet-spin and a small fraction with a triplet-spin 
ground state. Presence of a wide range of distributions of 
singlet-triplet energy gaps of dimerized V moments was 
claimed. The total gap distribution consists at least of 
three groups. Two of them, relatively narrow, are cen- 
tered at ~ 1 and 10 K and involve only 1/9 of vanadium 
moments. The third group, that corresponds to the 8/9 
part of V moments, is very broad and this leads to a 
wide range of spin couplings. Additionally at very low 
temperatures, Ri86 mK, a phase transition was observed 
in the temperature dependence of the ac-magnetic sus- 
ceptibility Xac(T). It was interpreted as a proximity to a 
quantum critical point at zero field. 
Several theoretical calculations of the electronic structure 
and exchange interactions in Na2V307 delivered conflict- 
ing results. A spin dimer analysis by Whangbo and Kooi 
lead to the conclusion that Na2V307 could be described 
by six mutually intersecting helical chains. A single he- 
lical chain has a spin gap and it is a good approxima- 
tion for magnetic properties of the whole system. This 
is not supported by experimental results. The ab initio 
Density Functional Theory (DFT) calculations by Saha- 
Dasgupta et ali^, used the N-th order muflin-tin-orbital 
(NMTO) based downfolding method to obtain a tight- 
binding Hamiltonian for V-3dxy orbitals. Their micro- 
scopically derived hopping integrals resulted in antiferro- 
magnetic (AF, J<0) intra-ring exchange integrals domi- 
nating the intcr-ring ones. This lead to the description of 
the system as weakly coupled nine-site rings with Ji=: - 
200 K and 32= -140 K intra-ring couplings. Mazurenko et 
al^ performed DFT calculations with the local density 
approximation (LDA)-t-J7 approach and suggested that 
due to the strong hybridization between filled and vacant 
3d vanadium orbitals the exchange interaction between 
V ions is dominated by strong ferromagnetic (F) inter- 
ring contributions, while the intra-ring exchange interac- 
tions are mainly antiferromagnetic. Their conclusion is 
that the resulting Heisenberg model is strongly frustrated 
and there is no simple way to predict the properties of 
the system. 

The main goal of our paper is to provide, by character- 
izing static and dynamic, nuclear and magnetic, micro- 
scopic and macroscopic responses, a solid experimental 
background necessary to derive a correct model of the 
magnetic exchange in Na2V307. 



II. EXPERIMENTAL DETAILS 

The Na2V307 polycrystalline samples used in this 
study have been prepared by a solid-state reaction 
method using Na4V207, V2O3 and V2O5 as described 
by Niitaka et al7 . The precursor Na4V2 07 was syn- 
thesised from a mixture of NaCOs and V2O5 with a 
stoichiometric ratio, on a gold sheet at 873 K under air, 
and kept at 473 K before use. V2O3 was reduced from 
V2O5 by 7% H2 + 93% Ar mixture gas. Under dry- 
and pure-Ar atmosphere, the stoichiometric amount of 



Na4V207-|-V203-|-V205 was well mixed, and pressed 
into pellets. The pellets were wrapped in gold foils, 
and loaded into a fused silica tube. The tube was 
subsequently sealed under vacuum and heated in a 
furnace at 978 K for a week. The product was checked 
as pure Na2V307 by x-ray powder diffraction^. 
Analysis of the crystal structure and its evolution with 
temperature was studied by x-ray single crystal and 
powder diffraction. A black fiber-like single crystal 
(10x10x700 /im) was extracted from the polycrystalline 
sample. X-ray diffraction studies have been performed 
at the BMOIA Swiss-Norwegian BeamHne (SNBL) of 
ESRF, France, at a wavelength of 0.722870 A, and at the 
Material Science Beamline (MSB) of SLS, Paul Scherrer 
Institute, Switzerland, with A=0. 620967 A. The crystal 
structure has been determined from a single crystal at 
SNBL at three different temperatures, 200 K, 80 K 
and 16 K, using an image-plate area detector MAR345 
and N2 or He flow cryostats to reach low temperatures. 
180 deg oj oscillation images with an increment of 1 
deg were collected with 1 min of exposure and the 
crystal-to-detector distance of 180 mm (see Table HI). 
The quality of the data at 16 K is lower than at 80 K 
and 200 K due to a lower resolution obtained with the 
complex low-temperature experimental setup. The data 
were corrected for the Lorentz-polarisation factor and for 
absorption. For the structure solution and refinements 
the programs SHELXS and SHELXL^ were used. 
Powder x-ray patterns were collected in a 0.5 mm quartz 
capillary at temperatures between 80 K - 470 K with 
the same setup used for the single crystal measurements. 
The experiment was performed on heating, with a heat- 
ing rate of 120 K/h. One exposure lasted 30 s and the 
capillary was rotated 30 deg/exposure for better powder 
average. For patterns recorded in the temperature 
range of 14 K - 130 K a Janis He cryostat at MSB and 
a microstrip detector were used with 30 s exposure. 
Structure refinements were performed using the Rietveld 
technique implemented in the FuUprof prograroiS. 
Magnetic excitations and lattice vibrations were in- 
vestigated on a 10 g polycrystalline sample. Inelastic 
neutron scattering measurements were performed on 
the neutron time-of-flight (TOF) spectrometer FOCUS 
at SINQ, PSI, Switzerland, at 1.5 K and 120 K. Three 
setups were exploited (Ai= 5.75, 4.85 and 1.7 A) to 
access energy-transfer ranges below 0.9, 2 and 20 meV, 
respectively, with optimal resolution. Contributions 
from an empty Al sample holder were measured and 
subtracted. 

A powder neutron diffraction experiment with full 
three-directional (XYZ) neutron polarisation analysis 
was performed on the D7 instruments^ at the Institute 
Laue-Langevin, France. The incident neutron wave- 
length 3.1 A was selected using a pyrolytic graphite 
monochromator. Polarisation of the incident and 
analysis of the scattered beam were performed using 
polarising mirrors. Helmholtz coils enabled rotation of 
the incident neutron spins along three mutually perpen- 
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dicular directions. The efficiency of the 132 detectors of 
D7 was determined using a standard vanadium sample, 
while the polarisation efficiency was determined using 
a quartz sample with entirely non-spin flip scattering. 
The measurements were performed at 1.9, 120 K (each 
60 h) and 300 K (2 h). Our experimental setup allowed 
us to measure the response of the sample within the 
wave vector range 0.5 < Q < 4 and to integrate 
the signal over the energy interval -co, 9 meV. The 
scattering from the sample was brought to the absolute 
scale using vanadium normalization. Absorption of the 
vanadium and the sample was taken into account. 
A powder non-polarized neutron diffraction experiment 
in the mK range was performed on DMC at SINQ, PSI, 
Switzerland, using a wavelength of A=2.457 A. The 
sample was mounted into the dilution insert of an ILL 
orange cryostat and two diffraction patterns, at 40 mK 
and 1 K, were collected. 

Measurements of the susceptibility on a polycrystalline 
sample were done at MPICPfS, Dresden in the 30 mK - 
3.5 K range in an impedance bridge using a SQUID as 
a null detector. In this arrangement, the amplitude of 
the ac field can be varied in fixed steps from 0.07 to 33 
mOe and the frequency can be chosen from 4 different 
values between 16 and 160 Hz. The sample was placed 
inside the mixing chamber of a dilution refrigerator. 
The residual field in the cryostat was less than 20 mOe. 
Specific heat measurements were performed using a 
quasi-adiabatic method in the temperature range 60 mK 
- 4 K at four values of apphed field 0, 0.25, 0.5 and 2 T. 
The decay of the remanent magnetization was measured 
at MPICPfS, Dresden. The sample was cooled in a field 
of 100 Oe from T= 200 mK to the desired temperature 
T<T/. Note that 200 mK is well above the critical tem- 
perature T/= 86 mK where a cusp in Xac occurs. After 
a waiting time of about 10 hours, the field was removed 
and the relaxation of the remanent magnetization was 
recorded with a fiux counter for about 20000 - 70000 
seconds. After these times, the sample was heated to 
200 mK and the expelled flux was recorded in order to 
obtain the total value of the remanent magnetization 

Mrem = Ma - M^. 



III. RESULTS 

A. Temperature evolution of the crystal structure 

In an attempt to understand the magnetic exchange 
in Na2V307 we performed detailed structural investiga- 
tions at low temperatures. The analysis of the 200 K 
x-ray single crystal diffraction data set suggested an 
apparent Laue symmetry 3m with mirror planes situated 
along the [100], [010] and [110] directions. No systematic 
absences of the space group P31c reported by Millet 
al} were observed. The direct methods were tentatively 
used with the following probable space groups P31m, 



FIG. 1: a) [001] projection of the Na2V307 structure. Ar- 
rangement of V atoms forming buckled rings (slices A, B) 
within the nanotube b) in the model of Millet et ali- and c) 
in our model. The three sites of vanadium are denoted by 
VI - red, V2 - pink, V3 - blue. Oxygen and sodium atoms 
are omitted for clarity, d) Four different families of exchange 
couplings Ji (grey), J2 (white), J3 (black), J4 (black dashed) . 
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P31m and P312, but no acceptable solution was found. 
Therefore, the presence of {010} twinning of the crystal 
with P3 or P3 symmetry has been considered. The 
centrosymmetric model presented in Table [Tl] with 50% 
of two twin domains explained well the apparent 31to 
symmetry. The R-values corresponding to the final 
refinements are presented in Table HI The anisotropic 
displacement parameters and occupancies of all atoms 
could be refined at 80 K and 200 K (see Table [ml). For 
the 16 K data set, due to its lower resolution, isotropic 
displacement parameters were used. 
Single crystal data collected at three temperatures 
indicate no structural phase transition. The structure 
resulting from our analysis is very similar to the model 
proposed by Millet et alX However, there are important 
differences. The arrangement of nine V atoms into the 
buckled ring (slice) B is the same as in the model of 
Millet (Figs. [T)d and c). But due to the absence of the 
glide plane in the P3 space group, the slice A' is oriented 
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FIG. 2: Thermal expansion of the lattice parameters a (red) 
and c (blue) from x-ray powder diffraction. Inset: The a/c 
ratio calculated from the SNBL (green), SLS (violete) powder 
patterns and SNBL (red) single crystal data. 
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differently compared to A. We identify the A'-B two-ring 
stacking unit as a structural building block of Na2V307. 
A detailed structural investigation resulted in precise 
values of the V-V distances and V-O-V angles presented 
in Tables IIVIVI Here a second index in the notation 
of V atoms is introduced to simplify the comparison of 
numerical and graphical presentations (see also Fig. [TJi). 
Based on these values it is possible to distinguish four 
different families of effective exchange paths. They 
are: 1** and 2"^ intra-ring and 1** and 2"'' inter-ring 
ones. In a simplified model for the magnetic exchange 
in Na2V307 each family may be represented by a 
single exchange couphng (Ji — J4, see Fig. [TJi). The 
1** intra-ring exchange path is characterized by the 
average shortest V-V distance of «2.985 A and the 
V-O-V angles of w99 deg. Thus it is the most im- 
portant coupling. The 2'"^ intra-ring path has rather 
large V-V distances (w3.63 A) and V-O-V angles of 
«140 deg. The 2"'' inter-ring path has smaller V-V 
distances (wS.SS A) and almost the same V-O-V angles 
of f«132 deg. Therefore, based solely on the atomic 
arrangement, the 2"^^ intra-ring and 2"^^ inter-ring paths 
might be equally important. The 1** inter-ring path 
will certainly result in a 'separate' coupling as the V-V 
distance is ~3.37 A and the V-O-V angles are wll9 deg. 
DFT calculations^ performed on the crystal structure 
proposed by Millet et alX identify the 1*** and 2"*^ 
intra-ring couplings as the most important ones. And, 
according to these calculations, the 2"'' intra-ring and 
2"'' inter-ring couplings are significantly different due to 
the orientation of the V orbitals. As small differences 
in the atomic arrangement between our model and the 
model of Millet, can significantly change the strength 



FIG. 3: Thermal evolution of the V-V distances calculated 
from single crystal x-ray diffraction at 200, 80 and 16 K. Top: 
1"* intra-ring distances, bottom: 2"'' intra-ring (green), 1"* 
(red) and 2"^* (blue) inter -ring distances. 
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and even the sign of the magnetic exchange^, we revised 
the electronic structure of Na2V307 performing new 
DFT calculations (see Section ITVT) . 

In the new model of crystal structure, within the 
vanadium-oxide nanotubes there are two inequivalent 
Na sites, 1(a) and 1(b). They are only partially occu- 
pied. This leads to a nonstoichiometric composition 
Na2-2;V307 with x=0.17. The requirement of neutrality 
of the crystal inquires that the missing positive charge 
is distributed on the neighboring atoms; either some 
amount of V^"*" (up to 5.7%) or intermediate valence 
states should be present. These minority states, if local- 
ized on V, could serve as defects introducing randomness 
into the V-V exchange couplings. 

The displacement parameters U of all atoms, including 
sodium, are small, almost isotropic (see Table IIIip and 
decrease with temperature lowering, as in conventional 
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materials. Therefore, the assumptions used to identify 
the low frequency phonon mode at 88 cm~^ as the Na 
rattling modei^ does not hold. Our Na^ mean square 
displacement amplitude < > is approximately six 
times smaller than assumed by Choi et al.^^ and there- 
fore the characteristic frequency would be « 174 cm^^, 
which is incompatible with the observed hypothesized 
rattling mode at 88 cm^^. 

One more important outcome of our x-ray diffraction 
study is the temperature evolution of the lattice pa- 
rameters a and c. Powder patterns reveal that the 
lattice compresses gradually on coohng (Fig. This 
compression is very anisotropic, as visualized in the 
a/c ratio in the inset of Fig. O Whilst above 300 K 
the expansion coefficients Aa/a and Ac/c lie in the 
range of 10^^/K, which is typical for oxides; several 
anomalous points in a(T) and c(T) are observed at lower 
temperatures. Near 300 K the decrease of c slows down, 
close to 100 K with decreasing temperature, a(T) starts 
to increase. Unfortunately we cannot supply microscopic 
details of all these changes, but we reckon that the 
changes near 100 K correlate with the anomaly observed 
in the ^^Na-NMR spin-lattice relaxation rate^ and the 
88 cm~^ mode detected in optical measurements^. 
The variation of the V-V distances and V-O-V angles 
obtained from single crystal refinements (Tables IIVIVI 
and Fig. [3]) allows us to conclude that the dominant 
tendency is a reduction of the V-V distances and V-O-V 
angles having the largest projection along the tubes and 
an increase of those having a substantial xy-component. 



B. Neutron scattering: polarized diffraction and 
time-of-flight inelastic scattering 

We performed neutron scattering experiments to ob- 
tain information on the static and dynamic magnetic 
correlation functions in Na2-2;V307. This system has 
significant incoherent scattering aincoh due to the pres- 
ence of vanadium and sodium, aincoh is comparable to 
the magnetic scattering of V (cri„coft = l-47 barn/str. f.u., 
fmag=0.437 barn/str. f.u.), additionally absorption is 
significant. Therefore, a polarised neutron diffraction ex- 
periment was carried out on D7 to separate magnetic 
scattering from nuclear coherent, nuclear spin incoherent 
and background contributions. The paramagnetic scat- 
tering was obtained by averaging the contributions ex- 
tracted from the spin-flip (SF) and non-spin-flip (NSF) 
partial cross sections 

" 'Jmag _ r)r " "x ^ " "v _ q " " z j 

dildu dfldu! dilduj dVldio 

d'^rr d'^a'^S^ d'^a^^^ d'^a'^^^ 

" '^mag ^ 2[ ^ y I 2 ] (1) 

dVLduj dfldio dQduj dfldio 

where z is perpendicular to the scattering plane. Fig. [3] 
(top) presents total, nuclear spin incoherent and mag- 
netic signals from the sample. Apparently the mag- 



FIG. 4: Top: total, nuclear spin incoherent and magnetic 
scattering at 120 K measured on D7 with Ai=3.1 A. Inset: 
The (Q, to) window accessible during the experiment is shown 
by the dashed area. Bottom: Absolute magnetic scattering 
at 1.8, 120 and 300 K. 




netic signal is tiny. The amount of magnetic scattering 
summed over the measured Q-range is only 13 % of the 
expected paramagnetic scattering of Na2-a;V3 07 formula 
unit (f.u.) in this Q-interval calculated as 

amag = ^QlNSiS + l)iro-/ff{Q), (2) 

where | results from the powder average, iV=3 is the 
number of spin S=l/2 V^^ ions per f.u., (ro7) = - 
0.54-10^^^ cm and f{Q) is the magnetic form factor of 
V^"*". The weakness of the signal caused long counting 
times to obtain reasonable statistics in the diffraction 
mode, i.e. integrating over all energy window, and it was 
impossible to measure the (Q, uj) distribution of the mag- 
netic scattering in TOF mode of D7. 
Fig. |4] (bottom) presents magnetic contributions at 1.8, 
120 and 300 K. The signals are almost equal and fea- 
tureless within statistical errors. This leads to important 
conclusions. Firstly, the majority of magnetic intensity 
remains out of the integrated energy and momentum 
window of our experiment shown as an inset in Fig. [4] 
(top). At 1.8 K dominantly elastic and neutron energy 
loss events are measured, while at 120 K additionally the 
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FIG. 5: Comparison of the magnetic scattering at 2 K mea- 
sured on D7 (red) and the calculated paramagnetic signal of 
V^^ (black solid) , both in absolute scale. Calculated low-lying 
excitations of an isolated nine-member ring (violet dashed) 
and of an isolated dimer (blue dotted) are given in arbitrary 
units. 
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FIG. 6: INS spectra of Naa-a^VsOr at 1.5 K (blue) and 120 
K (red) measured on FOCUS with Xi=1.7 A. The intensity 
has been integrated over the momentum transfer range 2.1 < 
Q < 5.7 A"\ Inset: Q-dependence at 120 K for 2.6 < Q < 
3.4 A"^ (green) and 4.1 < Q < 4.9 A"^ (magenta). 
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neutron energy gain side up to « 2fcT=20 meV is acces- 
sible. In both measurements the magnetic scattering is 
very weak and it only slightly increases from 1.8 to 120 K 
in the low Q range, which implies presence of high-energy 
states becoming populated at higher temperature. The 
total paramagnetic scattering is not recovered even at 
300 K (Figs. mis]), therefore, most of the magnetic inten- 
sity is located at even higher energy transfers. 
Secondly, observation of weak magnetic intensity with no 
oscillation in the Q-dependence disputes the presence of 
isolated clusters. Any isolated magnetic cluster would 
have a characteristic oscillation in the Q-dependence, 
which arises due to the geometric ^"^^^^-^ term in the 
neutron cross sectioi*i^: 



^^«,f(g)E7<,.,(|<S||Tj||S'>r) + 
^siniQAR) 



QAR 



< S||Tj||S' >< S'||Tj,||S > 



(3) 



where AR — |Rj — Rj/ 1 is the distance between the 
ions i and j of the cluster, Tj - are irreducible tensor op- 
erators. Fig. [5] presents the expected Q-dependences for 
an isolated nine-member ring and an isolated dimer cal- 
culated using programs of Weihe^^ . For the nine-member 
ring the Hamiltonian is -ff = JiS^^]^(Sj-Si+i)+ J2Sij (S^- 
Sj) (for i=2 j=4:, 9, for i=5 j=3, 7, for i=8 j=l, 6) and 
the two intra-ring antiferromagnetic exchange constants 
are Ji=-160 K and J2=-90 K, the best estimate from our 
present susceptibility data (see Section |V|. 



The geometric term is summed over the excitations be- 
tween the 16 lowest cluster levels which could contribute 
to the integrated energy window of the D7 experiment 
and Ai?'s correspond to the V-V intra-ring distances de- 
termined in this work. For the dimer the Hamiltonian 
is H = J(Si • S2) with the antiferromagnetic exchange 
constant J chosen such that the spin singlet-triplet exci- 
tation lies in the energy window of the experiment, AR 
is the average intra-ring distance. None of the cal- 
culated Q-dependences is consistent with the measured 
one. 

Yet, our observations are compatible with two following 
pictures: i) Na2_a;V307 is composed of complex objects, 
i.e. spin clusters with different sizes and exchange cou- 
plings, and the Q-dependence of each individual cluster 
is blurt out; ii) nine-member ring spin clusters are cou- 
pled into tubes, this leads to dispersion of excitations 
and to smearing of oscillations in the measured energy 
integrated S(Q, lo). An intermediate picture, nanotubes 
segmented into finite fragments with different length, also 
explains the observations. 

Inelastic neutron scattering collected at two tempera- 
tures 1.5 and 120 K (Fig. in the TOF mode on FO- 
CUS is dominated by incoherent scattering. No distinct 
magnetic excitations have been observed up to 20 meV in 
agreement with the D7 results. The weak mode observed 
near meV is most probably due to phonons since its in- 
tensity increases with Q and temperature (Fig. [S] inset). 
It is found close to the 88 cm~^ mode observed by optical 
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measurementsi^. Our observations, however, does not al- 
low us to specify which atoms (V or Na) are involved in 
this mode. 



C. Magnetic properties in the mK temperature 
range 

1. ac- susceptibility and specific heat 

We performed new experiments in the mK temper- 
ature range to characterize the previously found'^ very 
low-temperature state of Na2-2;V307. First of all, our ac- 



FIG. 7: Temperature and frequency dependence of real x ^nd 
imaginary x" components of ac-susceptibility in mK temper- 
ature range. 
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susceptibility measurements confirm the sharp anomaly 
in x(T) in the form of a "cusp". The temperature de- 
pendencies of both, the in-phase x' and out-of-phase x" 
components display peaks near 76 mK (Fig. [7]), which is 
slightly lower than the reported valued. Although x' has 
only a very weak frequency dependence, the magnitude 
and shape of x" (Fig. [7] inset) are strongly frequency de- 
pendent. In particular, the temperature where the max- 
imum of x' occurs decreases slightly with decreasing fre- 
quency. These observations hint for a slow dynamics of 
the V*^ moments reminiscent of spin-giass-likc processes 
with a freezing temperature Tf—76 mK. 
Although the anomaly in x(^) at Tf signals a change in 
the magnetic properties of Na2V307, the state below T/ 
is not long-range magnetically ordered. This is consistent 
with powder neutron diffraction. The patterns collected 
at 40 mK and 1 K were identical, no magnetic reflections 



occurred at low temperature. 

Fig. [5] inset displays the low-temperature part of the 
magnetic specific heat Cm obtained by subtracting the 
lattice contribution from the measured signal Cm{T) 
shows no clear features near Tf. This important result 
shows that the number of degrees of freedom involved in 
the spin-freezing phenomenon must be rather small com- 
pared to the number of ions. The second important 
feature is a broad Schottky-like feature near 0.4 K, which 
shifts to higher temperatures with applied field. Approx- 
imately 1/13 of the magnetic entropy is associated with 
this feature. It is worth to recall that there is also a 
broad feature at 5 K, where at most 1/9 of the magnetic 
entropy is released as found previously^. 
And finally, there is a strong upturn in Cm/T with de- 



FIG. 8: Magnetic part of the specific heat divided by temper- 
ature, Cm/T, versus T for different fields. Inset: Cm versus 
T below 1 K. 




T[K] 



creasing temperature below 0.15 K (Fig. [S]). This upturn 
is well represented by Cn/T = D/T^. D is of the or- 
der of 2x10"'' to 4x10-4 (J (mole f.u. K)-^) and Cn 
is only weakly field dependent. Paramagnetic impuri- 
ties (i.e. only weakly- or non- interacting magnetic ions) 
cannot be the origin of this feature because their contri- 
bution would be strongly suppressed by a field of 2 T, 
contrary to the experimental observations. The impu- 
rities, however, could come from a nuclear contribution 
to the specific heat due to i) nuclear quadrupolar energy 
levels and ii) nuclear Zeeman levels. 
First let's roughly estimate the contribution from nu- 
clear quadrupolar energy levels. Nuclei with spin / and 
quadrupolar moment eQ exposed to static electric-field 
gradients, where the largest component is denoted as eg, 
contribute to the specific heat (per mol of the considered 
nuclei)!*: Cat = Dq x T'^, with 



Dq ^ 1 (2/ + 2)(2/ + 3) / e'gg 
i? ~ 80 2/(2/ - 1) \ kB 



(4) 
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R is the universal gas constant, is the Boltzmann con- 
stant and e is the electron charge. The most promising 
candidates to provide a substantial contribution to the 
specific heat are the ^^V nuclei which havei^ / = 7/2, a 
large Q — 5.2 x 10"^^ m^ and low point symmetry. From 
Eq. Hwith £)q = 4 X 10""* (J mole f.u. K^i) it follows 
6(7 « 4 X 10^'^ V/m. This value is about two orders of 
magnitude too large compared with experimental NMR 
studies of other V compounds^S. Therefore, the nuclear 
quadrupolar levels seem to be improbable reason of our 
observed results. 

Second we consider nuclear Zeeman levels. Assume a 
fraction / of magnetic V moments, which are at/near 
freezing. The direct contribution from these moments 
was already discarded previously. However, an indirect 
contribution via the hyperfine field coupling to the nuclei 
should also be taken into account. The nuclei are influ- 
enced by the on-site moment via the core polarization. 
This would result in a local field of the order of 100 T 
(per jiB oi onsite V moment)^*'. The magnetic nuclear 
contribution to Cat is Cat = Dz/T^, with 



FIG. 9: Temperature and frequency dependence of the nor- 



Dz ^ (2/ + 2) fjNhm 
R ■' 61 \ kg 



(5) 



This takes into account that there are 3 V ions per 
formula unit and assumes that only a fraction / of V 
ions are magnetic, -/n = 7.045 x lO'^ (rad (sT)-iT)i^ 
is the nuclear gyromagnetic ratio of ^^V and Hi is the 
local magnetic field. With / = 7/2 and Hi = 100 T this 
leads to / = 0.001. We conclude that a tiny amount 
of "frozen" V moments acting on their nuclei via the 
hyperfine field would explain the observed C jv (T) at low 
temperatures. 



2. Decay of remanent magnetization 

The dynamics of the V moments on a long-time scale 
at temperatures below Tf was probed by measuring the 
decay of the remanent magnetization Mrem = Mq — M^o- 
We recall that in these experiments the sample is cooled 
to a given temperature T<Tj in a fixed external mag- 
netic field H =100 Oe. After a waiting time of about lOh, 
H is turned to zero at t=0 and the time evolution of the 
magnetization Mt is monitored, keeping T fixed. The de- 
cay of Mt — Moo does not follow a simple functional form 
and is extremely slow (Fig. [S]) indicative of a spin-glass 
like behaviour. Characteristic time scales are of the order 
of IC* sec, which we define as the time where Mt — Moo 
reaches a fraction 1/e of its initial value. As shown in 
the inset of Fig. [5] top, Mrem grows exponentially with 
decreasing temperature. 

For different spin-glass systems the magnetization de- 
cays following different functional forms. The functional 
form of Mt — Moo is determined by the distribution of 
barriers that the system should overcome to reach the 



malized remanent magnetization mr 



Top: the 



same cooling field of 100 Oe has been used and measurement 
performed at different temperatures. Inset: Mrem [arb. units] 
versus T. Bottom: different cooling fields of 100 and 10 Oe 
have been used and measurement performed at 40 mK. Inset: 
A versus T'^ . 



T" 



T" 



T" 



T" 



30 mK 
40 mK 
50 mK 




time [s] 



equilibrium. Our data are well characterized by the so- 
called "enhanced power law" given h^^^^: 



Mt ~ Af c 
Mo - Mo 



= exp{-A[ln(cji)]^} (6) 



This functional form is rather general, and it reduces, 
for instance, to a "logarithmic" type-decay for the 
particular case oi /3 ^ 1, A << 1 and t not too large. 
The underlying model assumes that as temperature 
decreases there is a spontaneous formation of clusters 
of correlated spins with non-zero total spin, which grow 
randomly and independently from each other leaving 
a "fluid" of individual spins in between. To reach 
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the thermal equihbrium a spin-reorientation of these 
clusters must occur, and barriers, which depend on the 
cluster volume V, must be overcome to reach thermal 
equilibrium. This yields a hierarchy of relaxation times 
and results in the enhanced power-law for the decay of 
the magnetization. 

The above form of the decay does not depend on the 
mechanism that drives the relaxation, which may be 
either quantum tunneling (QT) or thermal activation 
(TA). In addition the three parameters describing 
relaxation A, uj, (3 do not independently vary with tem- 
perature. For instance, the attempt frequency uj may be 
assumed T— independent. In case of quantum tunneling 
the parameters A and (3 are also T— independent, 
resulting in a temperature independent relaxation. This 
is clearly not the case in our experiments (Fig. [5]). 
In the case of thermal activation lo is T— independent, 
but A and (3 are correlated; namely, A (x (T/Tf)^ 
—. To a good approximation A versus T^^ deduced 
from our fits (inset of Fig. [H] bottom) is a straight line. 
Therefore, our data suggest a classical type of spin-glass 
phenomena driven by thermal activation, even at the 
lowest temperatures. 

Our observations are not consistent with the 'super 
paramagnetic' picture describing the Mni2^, Visr^, 
Fei3^ cluster systems. There, a well defined and 
unique magnetic barrier leads to a single relaxation time 
resulting in a single exponential decay. 



IV. AB INITIO DFT CALCULATIONS 

We have performed ab initio DFT calculations for the 
new crystal structure determined in this work (see Table 
II) and applied the NMTO-based downfolding method 
in order to obtain the effective vanadium- vanadium hop- 
ping matrix elements. We considered a crystal structure 
with 100% occupation of the Na positions. The analysis 
of the Na density of states (DOS) contribution near the 
Fermi level shows that the effect of Na is almost negli- 
gible concerning the V hybridizations, and therefore sig- 
nificant changes of the electronic structure are not to be 
expected by reducing the Na content in the calculations. 
Variations in the hopping integral determination between 
the crystal structure proposed by Millet et alX^ and the 
present one are only observed in (i) the F** and 2"'* intra- 
ring hopping parameters which are now almost identical, 
enforcing the frustration effects in this system, and (ii) 
the inter-ring hopping parameters are slightly larger than 
those obtained previously, but remain still weaker than 
the intra-ring hopping parameters. The analysis of the 
orbital orientation in the 2"'' intra-ring and 2"'^ inter- 
ring paths shows that, eventhough distances and angles 
are similar for both cases, the inter-ring hybridizations 
are oi dd — 6 type explaining the weakness of the hopping 
integral in contrast to the mixed dd — a dd — tt nature of 
the intra-ring orbital hybridization. Still, the exchange 



coupling constants between Vi and Vj , Jij are very sensi- 
tive to ferromagnetic J^j^ and antiferromagnetic J:^^^ 
contributions, Jy = Jf/' + J^j'''^' . Calculation of J^^^ 
due to the hybridization of Ydxy orbitals, shows that the 
antiferromagnetic intra-ring exchange interactions are a 
factor of four larger than the antiferromagnetic inter- 
ring exchange interactions. The ferromagnetic contribu- 
tions Jij^^ , which can be estimated following the Kugel- 
Khomskii model^i^ are, on the contrary, of the same 
order of magnitude for intra- and inter-ring couplings, so 
that in conclusion, the intra-ring interactions Jmtra are 
dominantly antiferromagnetic, while the inter-ring inter- 
actions Jinter havc important ferromagnetic components 
as pointed out also in Ref. i6|. 



V. SUMMARY 

We presented new experimental results which allow 
to advance the understanding of the puzzling properties 
of the nanotube system Na2-sV307. Revisited crystal 
structure belongs to the P3 space group, while the 
apparent 3m symmetry is due to twinning. The bucked 
ring of nine V atoms, is the same as in the model 
of Millet et alr^, but the orientation of rings in the 
nanotube is different. The A'-B two-ring stacking unit 
is the basic structural building block. 
We found off-stoichiometry in the sodium composition 
that might introduce randomness into the exchange 
couplings or presence of up to 5.7% of V^+ defects. 
The lattice compression with decreasing temperature 
is very anisotropic, with several anomalous points that 
correlate with previous NMR and optical measurements. 
With decreasing temperature the dominant tendency 
is a reduction of distances directed mostly along the 
tube and an increase of those having substantial xy- 
component. No experimental support of a structural 
transition, spin dimerization and Na rattling was found. 
In the neutron scattering experiments we observed 
a weak magnetic signal in the energy window of our 
experiment from -20 to 9 meV, implying that the 
majority of magnetic excitations are located above 20 
meV. This is consistent with specific heat data, which 
indicate that at most 1/9 of the total magnetic entropy 
is released below 20 K, and with the susceptibility 
data, which monitor one order of magnitude moment 
reduction below 100 K. The magnetic scattering does 
not have a pronounced Q-dependence expected for spin 
clusters and, therefore, does not support the models of 
dimer formation or presence of isolated nine-member 
rings. The data might be explained, however, by the 
existence of clusters with different size and couplings 
or/and significant inter-clustcr couplings smearing out 
oscillations in S{Q,J uj). We note, however, that these 
scenarios leave unexplained the prominent peak in the 
temperature dependence of the spin-lattice relaxation 
rate around 100 K'^. That peak was used as part of 
the evidence that spin dimerization may occur in that 
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temperature range. But our present findings show that 
it is improbable. 

Our present and previous experimental observations al- 
low us to suggest the following picture: The Na2-a;V3 07 
system consists of segmented S=l/2 spin nanotubes with 
strong AF intra-ring but also significant inter-ring cou- 
plings. We speculate that segmentation happens due to 
the defects occurring because of the deviations from ideal 
stoichiometry. 

Though the dominant feature of susceptibility, the up- 
turn at 100 K, can be well reproduced by an isolated nine- 
member ring with Ji= -160 K and J2= -90 K (Fig. [TO)) . 
other experimental observations on Na2-xV307, i.e. neu- 
tron magnetic scattering, object a picture of isolated 
nine-member rings and require the presence of inter-ring 
couplings or more complex spin objects. The susceptibil- 
ity is dominated by strong intra-cluster couplings and is 
not a sensitive measure of the weaker inter-cluster cou- 
plings. A similar situation happens in coupled systems 
of trianglea^S and tetrahedra^i^. 

ab initio calculations point to larger intra-ring than inter- 
ring hopping integrals. The evaluation of the exchange 
integrals suggest strong antiferromagnetic inter-ring in- 
teractions and moderate intra-ring interactions mostly of 
ferromagnetic nature. 

As there is a small magnetic moment at low T, we in- 
fer that Na2-a;V307 belongs to the class of odd-legged 
S=l/2 spin tubes with a gapless excitation spectrum 
and main spectral weight located at high T. As shown 
by Liischer et al2^ frustrated inter-ring couplings might 
lead to the absence of a spin gap, which otherwise should 
be present in odd-legged tubes. A single crystal inelastic 
neutron scattering experiment is needed to verify our hy- 
pothesis, to obtain the inter-ring couplings and to explain 
the low-temperature features at 0.4 and 5 K in specific 
heat. 

The basic question, whether the ground state is mag- 
netic or non-magnetic, still cannot be answered unam- 
biguously. The origin of the 76 mK anomaly and of the 
observed slow dynamics in the mK range is not imme- 
diately evident. Our data strongly suggest a classical 
"spin-freezing" type phenomenon. But, given the small 
amount of spins involved (of the order of 0.1%) it is dif- 
ficult to ascertain whether this behaviour is intrinsic to 
the ideal nanotube or it is induced by segmentation of 
nanotubes. 
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TABLE I: Details of single crystal x-ray diffraction data col- 
lection and structure refinement for Na2-a;V307 at three dif- 
ferent temperatures. 



SFB/TR49. We thank for the expert assistance by Dr. 
F. Gozzo, Swiss Light Source and Dr. L. Keller, SINQ, 
PSI. 



T (K) 


200 


80 


16 


a (A) 


10.8827(5) 


10.8747(3) 10.8748(1) 


c(A) 


9.5506(2) 


9.5382(1) 


9.5275(2) 


V (A«) 


979.57 


976.86 


975.77 


20max{deg) 


59.75 


59.81 


42.99 


Rint% 


0.0182 


0.0153 


0.0423 


Refl. collected 


10485 


1744 


689 


Refl. unique 


1744 


1744 


689 


restraints 








18 


refined parameters 


113 


113 


69 


R, wR (J > 2cr(J)) % 


1.24, 3.99 


1.23, 3.78 


1.91, 5.16 


GooF on F'^ 


1.23 


1.15 


1.17 


largest diff. peak max/min 


0.33/-0.32 


0.38/-0.33 


0.31/-0.27 


and hole/e-A"^ 
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TABLE II: Atomic coordinates xyz and isotropic displace- 
ment parameters Uiso(-IO^A^) in Na2-:rV307 at 200 K. The 
space group is P3, Z=6, twin proportion is 49.90(9) and 
50.10(9)%. The Nal and Na2 sites are occupied only 45.3(11) 
and 51.8(11) %. Atomic parameters at 80 K and 16 K can be 
obtained as supplementary materiaL 



atom 


site 


X 


y 


z 




VI 




0.31549(3) 


0.10972(3) 


0.40031(3) 


0.354(7) 


V2 


6(7 


0.31643(3) 


0.149623(3) 


0.75026(4) 


0.316(6) 


V3 


6ff 


0.31391(3) 


0.18514(3) 


0.10094(3) 


0.390(7) 


Nal 


la 











0.8(1) 


Na2 


Ih 








0.5 


1.2(1) 


Na3 


2d 


1 
1 


2 


0.0701(1) 


0.79(2) 


Nal 


2d 




0.4332(1) 


0.95(2) 


Na5 


6g 


0.5074(1) 


0.4805(1) 


0.2492(1) 


2.07(2) 


01 


eg 


0.4780(1) 


0.1491(2) 


0.4380(1) 


1.02(2) 


02 


Gg 


0.4889(1) 


0.2194(1) 


0.7517(2) 


1.02(3) 


03 


Gg 


0.4760(2) 


0.3083(2) 


0.0636(1) 


1.15(3) 


04 


Gg 


0.2475(1) 


0.9508(1) 


0.7091(1) 


0.59(2) 


05 


Gg 


0.2432(2) 


0.0664(1) 


0.9329(1) 


0.65(2) 


06 


Gg 


0.2526(1) 


0.1701(2) 


0.5668(1) 


0.67(2) 


07 


Gg 


0.2715(1) 


0.2987(1) 


0.7911(1) 


0.52(2) 



TABLE III: Anisotropic displacement parameters 
Uiji-IO'^A'^) in Naa-^VaOy at 200 K. 



atom 


Uii 


U22 


U33 


U23 


Ul3 


U12 


VI 


0.36(1) 


0.40(1) 


0.30(1) 


-0.044(9) 


-0.07(1) 


0.19(1) 


V2 


0.28(1) 


0.39(1) 


0.25(1) 


-0.023(9) 


-0.02(1) 


0.15(1) 


V3 


0.42(1) 


0.46(1) 


0.31(1) 


0.011(9) 


0.06(1) 


0.23(1) 


Nal 


0.7(2) 


0.7(2) 


1.2(2) 








0.34(8) 


Na2 


1.1(1) 


1.1(1) 


1.7(2) 








0.54(7) 


Na3 


0.79(3) 


0.79(3) 


0.80(6) 








0.39(2) 


Nal 


1.00(3) 


1.00(3) 


0.87(6) 








0.50(2) 


Na5 


1.16(4) 


1.56(4) 


2.51(4) 


-0.40(3) 


-0.52(3) 


-0.05(3) 


Ol 


0.69(6) 


1.39(7) 


1.08(6) 


-0.42(5) 


-0.25(5) 


0.58(5) 


02 


0.54(6) 


1.30(7) 


0.89(5) 


-0.12(5) 


-0.03(5) 


0.22(5) 


03 


0.69(6) 


1.31(7) 


1.11(6) 


-0.16(5) 


0.16(5) 


0.25(5) 


04 


0.80(5) 


0.41(5) 


0.34(5) 


-0.04(4) 


0.03(4) 


0.14(5) 


05 


0.99(6) 


0.48(6) 


0.38(6) 


-0.04(4) 


0.11(4) 


0.31(5) 


06 


1.23(6) 


0.77(6) 


0.33(6) 


-0.04(4) 


0.04(5) 


0.74(6) 


07 


0.65(5) 


0.67(6) 


0.38(5) 


-0.02(4) 


0.02(4) 


0.43(5) 
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TABLE IV: Shortest V-V distances (A) at 200 K, 80 K, 16 
K and difference Ad between the values at 16 K and 200 K. 
The notations of Figure [2}1 are used. 



V-V distance 


200 K 80 K 


16 K 


Ad 




1"* intra-ring 






11-31 


2.9770(4) 2.9648(4) 


2.955(1) 


-0.022(1) 


21-32 


2.9851(4) 2.9908(4) 


2.996(1) 


+0.011(1) 


11-21 


2.9917(5) 2.9976(4) 


3.001(1) 


+0.010(1) 




2"'^ intra-ring 






21-31 


3.6109(4) 3.6008(4) 


3.596(1) 


-0.015(1) 


11-22 


3.6588(5) 3.6635(4) 


3.663(1) 


+0.004(1) 




1^* inter-ring 






11-23 
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3.365(1) 
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24-31 
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2"'' inter-ring 






31-33 


3.5447(4) 3.5436(4) 


3.545(1) 


0.000(1) 


11-12 
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TABLE V: V-O-V angles (deg) at 200 K, 80 


K, 16 K as- 


sociated with the shortest V-V distances. The notations of 


Figure [2ji are 


used. 






V-O-V angle 


200 K 80 K 


16 K 




1°* intra-ring 






11-07-31 


98.01(6) 97.28(5) 


97.1(1) 


11-04-31 


99.56(6) 99.26(5) 


98.9(1) 


21-04-32 


99.13(6) 99.55(5) 


99.9(1) 


21-05-32 


101.31(6) 101.39(6) 


101.4(2) 


11-07-21 


98.19(5) 98. 


24(5) 


98.4(1) 


11-06-21 


101.86(6) 102.07(6) 


102.3(2) 




2"'' intra-ring 






21-07-31 


136.15(7) 135.41(7) 


135.2(2)) 


11-04-22 


142.50(7) 143.41(7) 


143.9(2) 




1°* inter-ring 






11-06-23 


119.16(7) 119.02(7) 


119.0(2) 


24-05-31 


119.56(7) 119.48(6) 


119.3(2) 




2""* inter-ring 






31-05-33 


132.20(7) 132.00(7) 


132.0(2) 


11-06-12 


133.49(7) 133.61(7) 


133.5(2) 
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